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ABSTRACT 

The 3D structure of an active region (AR) filament is studied using nonlinear 
force-free field (NLFFF) extrapolations based on simultaneous observations at a 
photospheric and a chromospheric height. To that end, we used the Si I 10827 A 
line and the He I 10830 A triplet obtained with the Tenerife Infrared Polarime- 
ter (TIP) at the VTT (Tenerife). The two extrapolations have been carried out 
independently from each other and their respective spatial domains overlap in 
a considerable height range. This opens up new possibilities for diagnostics in 
addition to the usual ones obtained through a single extrapolation from, typi- 
cally, a photospheric layer. Among those possibilities, this method allows the 
determination of an average formation height of the He I 10830 A signal of ~ 2 
Mm above the surface of the sun. It allows, as well, to cross-check the obtained 
3D magnetic structures in view of verifying a possible deviation from the force- 
free condition especially at the photosphere. The extrapolations yield a filament 
formed by a twisted flux rope whose axis is located at about 1.4 Mm above the 
solar surface. The twisted field lines make slightly more than one turn along the 
filament within our box, which results in 0.055 turns/Mm. The convex part of 
the field lines (as seen from the solar surface) constitute dips where the plasma 
can naturally be supported. The obtained 3D magnetic structure of the filament 
depends on the choice of the observed horizontal magnetic field as determined 
from the 180° solution of the azimuth. We derive a method to check for the 
correctness of the selected 180° ambiguity solution. 

Subject headings: Sun: filament — Sun: activity — Extrapolations. 

1. Introduction 

Active region filaments are observed above polarity inversion lines (PIL) where the 
vertical component of the magnetic field changes sign separating the two opposite polarities 
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(Babcock & Babcock 1955). They are called filaments when observed on the solar disc 



and prominences when observed above the solar limb, though the two terms refer to the 



same phenomenon (see Demoulin 1998; Mackay et al. 2010, for a review on filaments and 
prominences). Filaments/prominences are formed by plasma which has been lifted up above 
the solar surface and is cooler and denser than its surroundings. The necessary force to 
sustain this plasma is of magnetic origin. Here we shall focus on active region (AR) filaments, 
which are characterized by a strong horizontal component of the magnetic field along the 



filament axis reaching several hundred gauss (Kuckein et al. 2009, 2011 Guo et al. 2010 



Canou & Amari 2010 Jing et al. 2010). These filaments typically lie low above the solar 



surface (Lites 2005) compared to quiescent filaments. 



There are two classes of models that aim at describing the magnetic structure of fil- 
aments: the sheared magnetic arcade model, and the twisted flux rope (TFR) one. The 
former involves photospheric motions of magnetic footpoints near the PIL produced by var- 
ious mechanisms such as photospheric flows, vortical motions or the emergence of the upper 



part of a flux rope (Antiochos et al. 1994; DeVore & Antiochos 2000 Aulanier et al. 2002 



Welsch et al.]2005| |DeVore et al.|2005| >. The flows force the magnetic field to get reconfigured 



creating a variety of field configurations, as dipped field lines (Mackay et al. 2010). 



The twisted flux rope model assumes that the filament is formed by a dipped flux rope 



where the material can naturally get located and lifted up ( 


van Ballegooijen & Martens 


1989 


Leka et al. 


1996 


Titov &; Demoulin 


1999 


Amari et al. |1999; Lites 


2005 


Lites et al. 



2010). Flux emergence of a twisted flux rope to the solar atmosphere has also been nu- 



merically mo 


delled (Magara & Longcope 


2003; 


Fan & Gibson 


2004; 


Archontis et al. 2004 


Amari et al. 


2004 2005 Galsgaard et al. 


2005 


Cheung et al. 


2007 


Martmez-Sykora et al. 


2008 


Tortosa-Andreu & Moreno-Insertis 


2009 


Fan 2009 MacTaggart & Hood 2010). From 


observations, 


Okamoto et al. (2008, 2009) 


have shown that the change in the horizontal field 



from normal to inverse configuration can be interpreted as the signature of an emerging flux 
rope. This has been shown as well in 3D MHD simulations of an emerging flux rope, where 
the spectropolarimetric signal exhibits similar changes of the horizontal field from normal to 



inverse (Yelles Chaouche et al. 2009a) using the same pair of Iron lines Fe I 6301 A and Fe I 



6302 A as the ones of Hinode SP/SOT. 

The numerical models aiming to simulate the properties of filaments assume a certain 
number of conditions, like, e.g.: the mere pre-existence of a flux rope in the convection zone 
or in the atmosphere; the amount of twist necessary for a given flux rope to emerge into the 
atmosphere; its morphology and emergence velocity; etc. These assumptions and conditions 
can hardly be directly tested using observed magnetograms only, and it turns out to be very 
helpful to use extrapolations in order to get information on the observed 3D magnetic field 
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and its properties. This gives a feed-back to numerical models and allows fine tuning them. 
In recent years, extrapolations have played an important role in probing active region 



filaments properties (Guo et al. 2010 Canou & Amari 2010 Jing et al. 2010). This is due 



to (a) the progress in mathematical and computing techniques of extrapolations (e.g. Amari 
et al.||1997t |Wiegelmahnl|2004| |Wiegelmann et al.||20"06l |Schrijver et al.||2006| |Aly fc Amari 



2007 Metcalf et al. 2008) and (b) to the large improvement in spectropolarimetric facilities 
allowing to measure the full stokes vector with ever higher sensitivity and resolution (E.g. 



SP/SOT onboard Hinode, TIP at the VTT, THEMIS/MTR, etc) see e.g. Lopez Ariste et al. 



(2006); Okamoto et al. (2008); Kuckein et al. (2009); |Lites et al. (2010), and the review by 



Mackay et al. (2010) for further references. 



In the past, extrapolations of active region filaments have exclusively been computed 
using photospheric fields as boundary values. Knowing that the photospheric plasma is not 
completely force free, it is important to compare these extrapolations with ones of the same 
region but using chromospheric magnetic field as boundary values for the extrapolations. 
The plasma in the chromosphere is usually much smaller than in the photosphere, so that 
the magnetic field in the former must be closer to force-free than in the latter. Simultaneous 



observations of the magnetic field vector at the photosphere and chromosphere (Kuckein 



et 


al.||2011 




Sasso et al.||2011|) 


in 


Kuckein et al. 


( 


2011), (see 



see also Kuckein et al. 2009, 2010) we will perform extrapolations 



starting from the photosphere and, independently, from the chromosphere. This allows to 
cross check the extrapolation results, and test whether the latter have been affected by e.g. 
the finite plasma (3 and the preprocessing procedure. 

We will focus here on two representative times (one for each day of observation). Then 
for each of them we have two snapshots: one at the photosphere using the Si I 10827 A 
line and one at the chromosphere using the He I 10830 A triplet. The time evolution of 
the filament over the entire observation period can be further retrieved from the paper by 



Kuckein et al. (2011). 



2. Observations and data analysis 



In this paper we study a filament located along the polarity inversion line of active 
region NOAA 10781, which was found to be in its slow decay phase when the observations 
were performed. The full Stokes spectropolarimetric data were acquired with the Tenerife 



Infrared Polarimeter (TIP-II; Collados et al. 2007) at the German Vacuum Tower Telescope 



(VTT, Tenerife, Spain) on the 3rd and 5th of July 2005, at coordinates N16-E8 and N16W18, 
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respectively. Several scans were taken with the slit (Q/!5 wide and 35") parallel to the filament 
obtaining maps with a field of view (FOV) of ~ 26 x 25 Mm and ~ 22 x 25 Mm, for July 
3rd and 5th respectively. It is important to note that the FOV was not centered on the 
same location on both days. However, the upper half of the map from July 3rd overlaps 
with the lower half of the map from July 5th. The observed spectral range comprises the 
photospheric Si I 10827 A line and the chromospheric He I 10830 A triplet, which allows us 
to study simultaneously both heights and their magnetic coupling. The reader is referred to 



the papers of Kuckein et al. (2009, 2011) for an extensive description of the data and the 



magnetic evolution of this active region filament. 

Figure [T] shows slit-reconstructed maps of the two selected data sets. From left to right 
different wavelengths are presented: continuum, Silicon line core intensity, Helium red core 
intensity and Silicon Stokes- V magnetogram (at —150 mA from line center). The upper map 
was taken on the 3rd of July 2005, between 14:39-15:01, and clearly shows the filament in 
the He I absorption panel (third panel starting from top left). The filament extends along 
the vertical axis and lies on top of the polarity inversion line, as shown in the upper Stokes- 
panel. The lower map was observed on the 5th of July 2005, between 8:42-9:01. The spine 
or filament axis is still seen in the lower half of the He I panel, but above, a more diffuse 
filament is located on top of the pores and orphan penumbrae which have newly formed. 

Two different inversion codes were used to fit the Stokes profiles. A binning in the 
spectral and spatial domain was applied to increase the signal-to-noise ratio. Hence, the 
final spectral sampling is ~ 33.1 mApx -1 and the average pixel size is ~ 1". For the Helium 



10830 A triplet we used a Milne-Eddington inversion code (MELANIE; Socas-Navarro|2001 ) 
and for the Silicon 10827 A line we used the SIR code (Stokes Inversion based on Response 



functions; Ruiz Cobo & del Toro Iniesta 1992). Kuckein et al. (2011) showed that the 



non local thermodynamical equilibrium (NLTE) effects of the Si I 10827 A line reported by 



Bard & Carlsson (2008) do not significantly affect the inferred vector magnetic field from 
the SIR inversion code and therefore were not taken into account. We transformed the 
vector magnetic field from the line of sight (LOS) into the local solar reference frame. The 



180° ambiguity was solved using the AZAM code (Lites et al. 1995). Response functions 



to magnetic field perturbations for the inferred atmospheres were computed at different 
positions inside and around the filament. An average optical depth of logr ~ —2 for the 
formation height of Silicon has been obtained for both days. 
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Fig. 1. — Two snapshots: The upper row (3 rd of July); The lower row (5 th of July). From left 
to right: Continuum intensity, Si I 10827 A line core, He I 10830 A red core, and Stokes-V 
at -150 mA from Si I 10827 A line center. As an example, the location of the spine (i.e. the 
filament's axis) can be seen in the upper He I panel as an elongated darkening along the 
y-axis. This can also be seen in the lower He I and Si I panels. The dashed-line rectangles 
indicate the approximate overlap of the filament between the two snapshots. 
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3. Force-free field Extrapolation 



Using the photospheric or chromospheric magnetic field as boundary values, it is possible 
to calculate the 3D magnetic field vector in the atmosphere using force-free field extrapo- 
lations (FFF). The force- free (or zero-Lorentz force) assumption is equivalent to assuming 
that the electrical current j and the magnetic field B are parallel, or, using Ampere's law 
and cgs units: 



V x B 



Air , 



aB ; 



T 



a in Eq. [T]is the so-called force-free parameter and measures the level of twist of the field lines. 
Given that B is solenoidal, a must be constant along each field line. We use here nonlinear 
force-free extrapolations, where a is constant along each given field line but can change from 
one field line to another. The extrapolations are carried out using the optimization method 



(Wiegelmann 2004) 



The observed magnetic field is not necessarily force free. This is especially the case 
in the photosphere where the plasma (3 is not small (with /3 the ratio of gas to magnetic 
pressures). For instance 1 in the lower part of the photosphere (Yelles Chaouche et al. 



2009b, from calculations using 3D MHD simulations). 



Schrijver et al. 


( 


2006, 


2008 


): 



plemented a variety of tests to compute the coronal magnetic field using several NLFFF 
extrapolation codes. The different codes were tested using analytical, numerical, and obser- 
vational models. Among other conclusions, these papers have addressed the issue that the 
NLFFF extrapolations are sensitive to the boundary conditions (usually, photospheric vector 
magnetograms). These might suffer from uncertainties due to low signal-to-noise conditions, 
inaccuracies in the resolution of the 180° ambiguity, or more importantly from non-force-free 
conditions of the photospheric plasma. 

In order to make the observed magnetic field maps more consistent with the force-free 
assumption and greatly enhance the correctness of the force-free extrapolations it is necessary 



to apply some preprocessing (Wiegelmann et al. 2006 Metcalf et al. 2008). This consists 



of minimizing the magnetic force and torque (Aly 1989), minimizing the difference between 



the observed magnetic field and the preprocessed one, and applying a smoothing operator 
to remove the small scale variations of the magnetic field. The preprocessing routine used 



here is the one developed by Wiegelmann et al. (2006). It is useful to quantify the change 



produced on the magnetograms by the preprocessing procedure. For that we use the usual 



vector correlation (Schrijver et al. 2006) between the observed vector magnetic field and the 



preprocessed one. It is found that for the case of the Si I 10827 A vector magnetoeram the 
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correlation between the original and preprocessed field is approximately 0.84. In the case of 
the He I 10830 A vector magnetoeram, the correlation is about 0.90 between the observed 
and preprocessed magnetic field. This indicates that the observed He I 10830 A magnetic 
field is closer to a force-free condition than the Si I 10827 A case. 

The observed vector magnetograms have a relatively small field-of-view and are not 
isolated from the rest of the active region. This might introduce some errors in the NLFFF 
extrapolated model. Some authors proceed in enlarging the field-of-view by embedding the 
observed vector magnetograms in MDI maps, nevertheless this might introduce inconsisten- 
cies since the MDI magnetograms contain only the line-of-sight component of the magnetic 
field, and therefore would inconsistently connect magnetic field lines throughout the NLFFF 
model. It is therefore safer to just use the existing vector magnetograms (after preprocess- 
ing) as lower boundary for the NLFFF extrapolations. A further test of the effect of reduced 



field-of-view on the extrapolation results is carried out in Sec. 4.6 



The extrapolation starts with potential lateral and top boundaries, whereas the bottom 
one is the preprocessed observed vector magnetogram. A description of the code can be found 



in (Wiegelmann 2004 Schrijver et al. 2006 Metcalf et al. 2008). The NLFFF extrapolated 



model is expected to be more correct in the lower central part of the computational domain 



(Schrijver et al. 2006) 



The extrapolations are computed in a box of 33 x 38 x 38 grid points in the x, y and 
z directions respectively, x and y being the horizontal directions perpendicular and parallel 
to the filament's axis respectively. This yields a distance of ~ 650 km between grid points. 
The same distance between grid points is used in the vertical direction. 



4. Results 

4.1. Extrapolations starting from the photosphere 

The vertical component of the magnetic field (Bjj h ) inverted from the photospheric 
Si I 10827 A data is plotted in Figure [2] using a color scale. The lower case "z" in Bf l '° 
refers to the vertical component of the magnetic field, while the upper case "O" and "Si" 
refer to the field obtained from Observations using the Si I 10827 A line. The arrows in 
Fig. [2] indicate the direction of the horizontal magnetic field. The length of the arrows is 
proportional to the horizontal field strength. In the vicinity of the polarity inversion line, 
some arrows point from the negative polarity to the positive one, which is called inverse 
configuration. A normal configuration is obtained when the field arrows point from the 
positive to the negative polarity. The inverse configuration in the PIL region suggests that 
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x [Mm] 



Fig. 2. — Background : Vertical component of the magnetic field on July 5 at the photo- 
sphere (-Bf -°). The arrows represent the horizontal component of the magnetic field. The 
field strength is proportional to the length of the arrows. 



Fig. 3. — Top: 3D view of the magnetic field of the filament on July 5 . The extrapolation 
results are computed using the observed photospheric magnetic field as lower boundary 
condition. Top background image: vertical component of the magnetic field at the bottom 
of the extrapolation cube. Purple and blue: B^ hE > 0; green and red: B^ l,E < Bottom: 
Extrapolation results using the chromospheric magnetic field on July 5 th as lower boundary 
condition. Bottom background image: the vertical component of the magnetic field at the 
bottom of the extrapolation cube B^ e ' E (same color scale as the top panel). The field lines 
are plotted in the filament region. The field lines that cross the lateral sides of the box are 
cut 2 pixels away from the boundaries to avoid plotting the part that might be affected by 
the lateral boundaries. 
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the filament above is formed by a twisted flux rope. 



The extrapolated magnetic field is plotted in the top panel of Figure [3] using UCAR's 
Vapor 3D visualization package (www.vapor.ucar.edu). This 3D view shows a sample of 
magnetic field lines in the filament region in order to outline its structure. The background 
map represents the vertical component of the magnetic field {Bf l)E ) at the bottom of the 
extrapolated cube. The upper case "E" and "Si" refer to the field obtained by extrapolations 
using the Si I 10827 A line vector magnetograms as lower boundary condition. In the 
background map of the top panel of Figure |3j purple and blue colors correspond to B® %)E > 
0; green and red are reserved for B^ t,E < 0. The magnetic field forming the filament is 
distributed as a flux rope with dips in the lower half of the figure (i.e., with field line 
stretches which are convex as seen from the solar surface). This goes in the same direction 



as the findings of Guo et al. (2010); Canou & Amari (2010) and Jing et al. (2010). 



Bard & Carlsson (2008), have studied the properties of the Si I 10827 A line in NLTE 



conditions. They have determined using empirical models that the average formation height 
of Si I 10827 A is about 320 km for a sunspot umbra and about 541 km in the quiet sun. In 
the case of the observed filament, we expect the average height of formation to be between 
these two boundary values. 

At this height of formation the photospheric plasma does not have small values of the 



plasma (5. Nevertheles, there are theoretical and observational indications (Wiegelmann 



et al. 2010 Martinez Gonzalez et al. 2010) that force-free extrapolations lead to a good 



retrieval of the actual magnetic field at higher atmospheric layers. 



4.2. Extrapolations starting from the chromosphere 

We have seen in the previous section, that the extrapolation studies have been conducted 
exclusively using photospheric spectral lines. Here we would like to present an extrapolation 
using the Helium triplet He I 10830 A as lower boundary for the extrapolations. Chromo- 
spheric lines as H a and He I 10830 A have been used to identify filaments but never as 
boundary condition for extrapolations. 

The bottom panel of Figure [3] shows a 3D display of a sample of field lines at the filament 
region. In this case the field lines harbor some shear in the lower half of the image but with no 
twist. Field lines have a normal configuration (from positive to negative), except at the very 
central part of the filament (shown by the yellow field line) where the field lines are almost 
parallel to the filament axis. The bottom panel of Figure [3] suggests that the magnetic flux 
rope forming the filament (Lower panel of Fig. p| is located below the height of formation 
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of the He I 10830 A triplet. We will come back in detail to this point in the next sections. 

4.3. Height of coincidence between the photospheric and chromospheric 

magnetic fields 

A natural question to ask when analyzing the extrapolations starting from the photo- 
sphere and the chromosphere is: does the magnetic field from the two extrapolations agree 
with each other? In order to answer this question, we first plot the field lines that result 
from the two 3D extrapolation boxes (Fig. [1]). The field lines start inside a 2 Mm deep box 
which extends along most of the filament. The field lines can develop outside of the box but 
will be cut if they go below it. This configuration allows selecting field lines above a suitable 
height in order to make it possible to compare field lines starting from the photosphere with 
those starting from the chromosphere. 

Therefore, the field lines in the left panel of Fig. [4] are all above an imaginary boundary 
situated at w 1.6 Mm above the bottom of the extrapolation cube. We have tried various 
heights and found that this one (~ 1.6 Mm) gives a good agreement between the field lines 
from the photospheric extrapolations (left panel) and the ones from the chromospheric ex- 
trapolations (right panel of Fig. [4]) which are drawn directly from the lower chromospheric 
boundary. A visual inspection of both panels suggests a good agreement between the mor- 
phology of the magnetic field obtained from photospheric and chromospheric extrapolations. 

To provide a quantitative comparison between the two sets of extrapolations, we com- 
pute the magnetic field strength in the photospheric extrapolation cube and look for the 
height where this field strength matches the observed chromospheric one. This operation is 
repeated along each column of the photospheric extrapolation cube. At the height (H) where 
the two field strengths coincide, the magnetic field components B^ l,c , By t,c , B^ h are stored 
as 2D arrays (The uppercase "C" stands for Coincide). For instance, the horizontal magnetic 
field from the photospheric extrapolation (B^ t,c , By l,c ) across the field of view is plotted in 
Fig. [5] and represented by red arrows. The one from the observed Helium is plotted as black 
arrows. The arrows length is proportional to the horizontal field strength. The top panel of 
Fig. [5] shows a relatively good agreement between the two magnetic fields indicating that the 
observed magnetic field B He, ° is well matched by the extrapolated B Sl)E field especially at 
the filament region. This goes in favor of validating the Silicon photospheric extrapolations. 
The average height where B St,E matches B He, ° is 1.57 Mm in the filament region, with a 
standard deviation of 0.66 Mm. Therefore, it is possible to determine an average formation 
height of the He I 10830 A triplet. This would be 1.57 Mm + the formation height of the Si I 
10827 A which leads to approximatively 2 Mm above the solar surface (r = 1). 
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Fig. 4. — Left panel: Field lines from the photospheric extrapolation. They are drawn above 
a minimum height of ~ 1.6 Mm above the photosphere (to make them comparable with 
the chromospheric magnetic field). Right panel: Magnetic field lines from the chromospheric 
extrapolations. They are drawn directly from the bottom of the corresponding extrapolation 
box (already at chromospheric height). The two sets of field lines start inside a 2 Mm deep 
box placed along the filament. Field lines are allowed to end outside of the starting box, but 
are cut if they go below it. The background images are as described in Fig. |3j and the data 
are from July 5 th . 
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Fig. 5. — Red arrows: Flow representation of the horizontal magnetic field from the pho- 
tospheric extrapolations at the height where the extrapolated photospheric field strength 
matches the observed field strength in the Helium vector magnetogram. Black arrows: Flow 
representation of the observed horizontal magnetic field in the He I 10830 A vector mag- 
netogram. Top panel: Solution obtained with the correct 180 degree ambiguity solution 
computed with the AZAM utility (Lites 2005, and references therein). Note that the solu- 
tions are obtained for both photospheric and chromospheric fields. Bottom panel: Similar 
to the top panel but using the other solution of the 180 degree ambiguity. All data sets were 
taken on July 5 th . 
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4.4. A complementary method for solving the 180 degrees ambiguity 

The top panel of Figure [5] indicates that the magnetic field extrapolated from the pho- 
tosphere matches well the one observed at the chromosphere. It is nevertheless interesting to 
redo the calculations but for the case where the observed horizontal magnetic field is taken 
from the solution with 180 degrees shift of the azimuth in the LOS frame. This is done to 
test the coherence of the results with the other magnetic field configuration. 

From the lower panel of Figure [5] it is found that, in this case, the magnetic field 
extrapolated from the photosphere and the one observed at the chromosphere have a quite 
different orientation in the horizontal plane although they have similar field strengths at the 
height of coincidence. This clearly discards the solution with 180 degrees shift and provides 
a complementary method to test the validity of the chosen horizontal field as a result of the 
various routines used to solve the 180 degrees ambiguity of the horizontal magnetic field. 



4.5. Location of the dips along the filament field lines 

The material carried by a filament must be sustained against gravity by the magnetic 
curvature force associated with the dips in the field lines. Using the extrapolated magnetic 
field from the Si vector magnetogram, we compute all the field lines that have a dip of at 
least 650 km depth (equivalent to the size of one grid cell) . Figure [6] indicates that such field 
lines with dips exhibit a twisted structure, revealing the topology of the filament's magnetic 
field. This X-Z view shows that the whole structure is lying relatively low. The axis of the 
filament is located some 1 Mm above the formation height of the Silicon line. This indicates 
that the axis of the TFR forming the filament is located below the formation height of the 
Helium triplet. The lower part of the field lines in Figure [6] exhibits an inverse configuration 
(from negative to positive) whereas the upper part has a normal configuration. Recall that 



the average formation height of the Helium triplet (Sec. 4.3) in the filament region is 1.57 
Mm above the formation height of the Silicon line. At that height the field lines have a 
normal configuration. This helps understanding the lower panel of Fig. [3j where the Helium 
extrapolations exhibit no twist, which also fits with the presence of a normal configuration 
magnetic field at the lower boundary of the extrapolation box. It can be seen from Fig. [6]that 
the field lines make slightly more than one turn. These field lines extend over a horizontal 
distance of ~ 18 Mm. This leads to ~ 0.055 turns per Mm. 

Following each field line in Fig |6j it is possible to calculate the location of the lowest 
part of the dips. These places are likely to be the preferred location where the plasma 



inside the filament would be gathered. Guo et al. (2010) and Canou & Amari (2010) have 
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Fig. 6. — X-Z view of magnetic field lines having a dip of at least 650 km depth. The compu- 
tations use the 3D magnetic field extrapolated from the Si I 10827 A vector magnetogram. 
This reveals the morphology of the flux rope with field lines making more than a full turn. 
In order to understand the morphology of this twisted flux rope, it is useful to have a look 
at the top panel of Fig. [3] showing another perspective of the filament magnetic field lines. 
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found a correspondence between the location of dips and the filament as observed in H Q . 
It is reasonable to argue that in order to build up opacity (e.g. in the H a spectral region), 
filaments are filled with material which will be preferentially located in dips. 

Figure [7] shows the locations of dips on top of the Si I 10827 A line core intensity image 
(left panel) and the intensity image of the red core of the He I 10830 A triplet (right panel). 
There is an agreement between the location of the filament and the dips, which are located 
in the lower half of the image where the filament has a flux rope structure. There are no 
dips in the upper part of the panel since the magnetic field in this region has a normal 
configuration with no flux rope. An even clearer case is seen in Fig. [8] where we repeated 
the same calculations as for Fig. [7] but for a snapshot taken on the 3 rd of July. Note the 
correspondence with both line cores. This analysis also provides a further indication on 
the location of the filament (as the H a spectral signature of a filament), and allows also to 
corroborate the extrapolation results since the spectral signature of the line cores of Si I 
10827 A and He I 10830 A are independent of the extrapolation analysis which involves the 
magnetic field and the force-free condition. 



4.6. Testing the effect of the reduced field of view on the extrapolations 

The observations made with TIP-II have a somewhat reduced field of view. Ideally the 
field of view should include the whole filament or even the whole active region (in general, 
a flux balanced region). The fact that we use part of the filament might influence the large 
scale magnetic connectivity through the filament. It is very useful to make quantitative tests 



to clarify the effects of reduced field-of-view (private communication, Aulanier 2011). We 
present here a test in which we cut part of the observed filament and analyze how that affects 
the extrapolations. The chosen reduced field-of-view region is the one delimited by the two 
dashed lines in Fig. [2j The resulting region is the one shown in the top-right panel of Fig. [9] 
corresponding to the 5 th of July. Similarly, the reduced field-of-view portion for the 3 rd of 
July is shown in the lower-right panel of Fig. |9j On these two panels we have drawn the 
location of dips using yellow "U" shapes. For the left column of panels in Fig. [9] we used the 
extrapolation of the whole field-of-view of Si I 10827 A vector magnetogram, but the field 
line calculation to find the locations of dips start from the lower dashed line until the upper 
dashed line. This covers a region similar to the reduced field-of-view cases on both the 3 rd 
(bottom-left) and 5 th (upper-left) of July. Figure [9] exhibit less locations of dips than the 
cases of Figs. [7] and |HJ This is primarily due to the fact that the field lines in the reduced 
field-of-view portion are shorter than the ones in the original full field-of-view case, and 
therefore there are less field lines that fulfill the criteria of having dips of 650 km depth. A 




Mm Mm 

Fig. 7. — Location of dips (drawn in yellow "U" shapes). These dips are calculated from 
the photospheric extrapolations for both panels. Left panel: locations of dips on top of Si I 
10827 A line core intensity. Right panel: locations of dips on top of the red core of He I 
10830 A. The data on both panels are from the 5 th of July. 
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comparison between the left and right panels in Fig. [9] shows that they have some similarities 
indicating that decreasing the field of view alters only slightly the structure of the obtained 
filament and more precisely the location of dips. Figure 10 provides a closer look to the field 
lines forming the dips plotted in the upper panels of Fig. M The field lines in the upper/lower 



panels of Fig. 10 correspond to the upper left/right panels of Fig. [9] respectively. The two 
sets of field lines are very similar indicating that the fact of reducing the field of view does 
not severely affect the magnetic structure of the extrapolated filament. Nevertheless, a more 
complete test of the effect of field of view on the results of extrapolation of active region 
filament would be useful and interesting. That would include the observation of a whole 
filament (ideally in an isolated and flux balanced region) and repeated extrapolations while 
reducing gradually the field of view. In the present work we were limited with the original 
field of view, therefore we could not run such a test. We intend to include a test of that sort 
as part of a future investigation. This test is complementary to the more complete series of 
tests presented in (e.g. |Schrijver etaLl|2006[ |2008[ |Metcalf et aT]|2008| |De Rosa et aT]|2009[). 



4.7. Magnetic field increasing from the photosphere upwards 



Due to magnetohydrostatic equilibrium, the field of concentrated vertical magnetic 
structures usually decreases with height above the photosphere. This is a typical case in 
a stratified atmosphere. Nevertheless, vertical gradients of the field strength, indicating 
stronger fields in the upper layers, have been observed (e.g., in polar crown prominences as 
shown by Leroy et aL] [l983) and modelled in the past ( Aulanier fc Demoulin||2003 ) and have 
been ascribed to the presence of dips ( Anzer|[l969 ; Demoulin fc Prie"st||1989 ). Here we would 
like to analyze the possibility of an upward increase of the magnetic field strength in the 
filament region. The upper panel of Figure [IT] shows the difference between the observed 
field strengths of and \B\ He, ° . For the sake of clarity, the figure only shows negative 

values of those quantities. A proxy of the location of the neutral line is indicated by the 
contour plots of B^ e, ° (i.e. small values of B^ e '° coincide with the neutral line location). 
In the central part of the figure, the values of \B\ He, ° are larger than \B\ Sl '°, indicating 
that the magnetic field in the chromosphere is larger than the photospheric one, suggesting 
that the magnetic field has increased upwards. In order to understand this phenomenon we 
look for the height where the extrapolated magnetic field from the photosphere reaches its 



maximum along each column. This height is plotted in the lower panel of Figure 11 It can 
be seen in the central region of the figure that the magnetic field reaches a maximum well 
above the photosphere, at an average height of about 1 Mm above the formation height of 
the Si I 10827 A line. 
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Figure [12] schematically explains the situation observed in those pixels where the mag- 
netic field increases upwards. First the magnetic field increases up to a maximum reached 



at the height HO (plotted, in fact, in the lower panel of Figure 11), then the magnetic field 
decreases gradually until reaching the height of formation of the Helium signal (HI) where 
the magnetic field is still larger than the one at the photosphere. It is also seen that there 
are two solutions for the height of formation of the Helium signal (HI and H2). The solution 
H2 is ruled out because it leads to abnormally low formation height. The HI solution is kept 



and is consistent with the average height of formation of the Helium signal (see Sec. 4.3). 



The reason why the magnetic field increases lies in the fact that the corresponding part 
of the filament resembles the structure of a force-free flux rope. In an ideal force-free flux 
rope, the azimuthal component of the field yields an inward-pointing magnetic tension force 
which is balanced by an increase of the magnetic pressure toward the rope axis. This leads 
to an increase of the magnetic field strength in the central regions of the rope. 



4.8. Complementary comments on the filament observed on the 3 of July 

Along the previous sections we have used observations taken on the 3 rd and 5 th of July, 
but with some focus on the filament as observed on the 5 th of July. This has been done in 
order to convey a clear message and not overwhelm the reader with information about the 



time evolution of the filament (This time evolution will be discussed in Kuckein et al. (2011 )). 
Additionally, the filament on the 5 th of July is very low lying (the filament's axis is located 
below the formation height of He I 10830 A). This represents a peculiar and interesting case 
which we have chosen to address. In order to complete the picture of the actually observed 
filament, we will briefly describe some further features of the filament as observed on the 3rd 
of July. 

It can be seen in Fig. [8j that the line core of the Si I 10827 A and the red core of the He I 
10830 A show a spectral signature reminiscent of the axis of the filament on July the 3 rd . The 
locations of dips along the filament coincide well with the line core signature (Fig. [8]). The 
filament in this case is located higher in the atmosphere and therefore imprints a clear spectral 
signature on both photospheric and chromospheric lines. The observed magnetic field from 
photospheric vector magnetograms exhibits an inverse configuration (field lines going from 
negative to positive polarities) in the filament region. The observed chromospheric vector 
magnetograms indicate that the magnetic field in the filament region has a slightly inverse 
configuration with most of the field being parallel to the filament's axis. From extrapolations 



(Fig. 13), it can be seen that the filament is formed by a twisted flux rope. This is revealed 



from both photospheric and chromospheric extrapolations. 
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5. Discussion and summary 

This paper is dedicated to building diagnostics about an active region (AR) filament 
using non-linear force-free extrapolations starting from the photosphere and the chromo- 
sphere. Extrapolations of AR filaments are usually done using the photospheric magnetic 
field as a lower boundary condition. This is generally convenient since the photospheric 
magnetic field is more accessible to observational techniques for various reasons, e.g. i) it 
has a large magnetic field implying a large Zeeman splitting, ii) the photon flux coming 
from the photosphere is large producing a high signal to noise ratio which makes it easily 
detectable, iii) there are many available instruments measuring photospheric field, etc. It is 
nevertheless interesting to probe the properties of filaments from two layers (or multi-layer 
if possible) since this approach provides several advantages like testing the validity of the 
extrapolation solutions. It allows as well testing the consistency between the extrapolations 
starting from the photosphere and the chromosphere. We just began to see the potential 
of multi-layer extrapolations, e.g. the determination of the matching height between the 
extrapolated photospheric field and the observed chromospheric one. This can be used as a 
method to determine the relative height of formation of higher- lying lines/multiplets com- 
pared to lower-lying ones. We find that the average formation height of the He I 10830 A 
triplet is about 2 Mm above the surface of the sun. Multi-layer extrapolations provides also 
a complementary method for testing the 180° ambiguity solutions. 

The AR filament studied here, harbors a TFR structure with a rather low-lying axis on 
the 5 th of July, where the axis is located at about 1.4 Mm above the solar surface. This height 
is below the formation height of the Helium signal. This result is consistent with the fact 
that the magnetic field as observed in Helium magnetograms harbors a normal configuration. 

In the present study we have essentially focused on a snapshot taken on the 5th of July 
(see Figs. 1 and 2). We have also briefly discussed about a snapshot observed on the 3rd 
of July (Figs. 1 and 12). Recall that on the 3 rd of July, the filament's axis was higher 
located and clearly seen in the spectral signature of Si I 10827 A and He I 10830 A line 
cores. Actually, as discussed in section 2 (observations and data analysis), the observations 
used here belong to a time series of snapshots covering the 3 rd and 5 th of July 2005. In 
order to build a clear and detailed diagnostic we have essentially focussed on studying one 
case (on the 5th of July). A more detailed analysis about the temporal evolution of the 



filament will be addressed at a later stage (see also Kuckein et al. 2011). Nevertheless, it is 
worth commenting on the other snapshots. An essential thing to mention, is that snapshots 
observed on the same day exhibit similar properties (height of location of the filaments axis, 
magnetic field strength, etc). This means that the rest of the snapshots observed on the 
5 th of July, have quite similar properties to the one described here, and similarly for the 
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snapshots taken on the 3 of July. 



For the filament on the 5 th of July, we have seen that, if taken separately, the photo- 
spheric and chromospheric extrapolations would lead to conclude that on the one hand the 
filament is formed by a TFR (from photospheric extrapolations), whereas on the other hand 
there is no TFR (from the chromospheric extrapolations). The multi-layer extrapolations 
helped solving this apparent discrepancy by showing that depending on the location of the 
TFR it might be seen at one layer or the other while the two extrapolations remain coherent 
with each other. Let us assume hypothetically that an active region filament is such that 
the lower part of the TFR is located at the upper chromosphere and does not extend to 
the photosphere. In this case, it would be seen in the chromospheric extrapolations, but 
the photospheric ones would not necessarily reveal it, especially since the photosphere is a 
relatively high plasma-/? medium and therefore the surrounding magnetic field in the AR 
might considerably alter the magnetic field near the PIL making the horizontal component 
of the magnetic field harbor a normal configuration. If we would only observe the photo- 
spheric field, this would lead to finding no TFR in the filament region, whereas in fact the 
filament's structure is a TFR. In fact this is true only if the electric currents present in the 



region are not reminiscent of a TFR, since it has been shown by Valori et al. (2010) that 



NLFFF extrapolations have the potential to find a TFR even if its dips are not present at 
the magnetogram level. Nevertheless, this hypothetical scenario highlights the usefulness of 
multi-layer extrapolations to help probing the properties of filaments and extracting further 
information about the observed magnetic structure. 

In the case studied here, the material in the dense part of the filament, along its axis, 
is confined inside a flux rope where the field lines are bent producing an excess of magnetic 
tension towards the center of the TFR. This tension is compensated by the outward oriented 
magnetic pressure producing an enhancement of the magnetic field inside the TFR. This 
is reflected by a vertical positive gradient of the magnetic field from the photosphere up 
to about 1 Mm above the formation height of the Si I 10827 A signal. In brief, the field 
increases from the photosphere upward in the flux rope region. 

There are only few studies dealing with extrapolations of AR filaments, the most recent 



being from Guo et al. (2010); Canou & Amari (2010) and Jing et al. (2010). These studies 



picture AR filaments as TFR with the exception of Guo et al. (2010) who have shown that 



TFR and dipped sheared arcades can co-exist in the same filament. We also find a flux 
rope from photospheric extrapolations, but have seen that the chromospheric extrapolations 
might show a different picture if the observed spectral line/multiplet is formed above the 
filament's axis. There is a clear necessity of further observational campaigns with the ap- 
propriate extrapolations to further clarify the properties of AR filaments and their temporal 
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evolution. Multi-layer extrapolations would be helpful to compare the resulting 3D structure 



of filaments with the variety of theoretical models existing in the literature (Mackay et al. 
20101). 
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Fig. 9. — Location of dips calculated from the photospheric extrapolations. Top row: loca- 
tions of dips on top of Si 1 10827 A line core intensity on the 5 th of July. Lower row: locations 
of dips on top of Si I 10827 A line core intensity on the 3 rd of July. Left column: the extrap- 
olation of the whole field-of-view of Si I 10827 A vector magnetogram is used, but the field 
line calculation to find the locations of dips start from the lower dashed line until the upper 
dashed line. This covers a portion similar to the reduced field-of-view test. Right column: 
the extrapolations of the reduced field-of-view of Si I 10827 A vector magnetogram are used 
(same size as displayed in the background Silicon line core image). Dips are calculated in 
this reduced field-of-view. 
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Fig. 10. — Lower panel: field lines with dips calculated in the extrapolation with reduced 
field of view corresponding to the region indicated in the upper-right panel of Fig. [9j Upper 
panel: field lines forming dips computed from the extrapolation with the whole field of view, 
but calculated inside the region includes between the two dashed lines in the upper-left panel 
of Fig. [9j 
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Fig. 11. — Upper panel: Difference between the observed |_B| S *'° and \B\ > . The positive 
values are clipped to in order to clearly see the negative values. Lower panel: The height 
where the magnetic field strength from the photospheric extrapolations reaches a maximum. 
The data on both panels are from the 5 th of July. The location of the neutral line is 
approximately indicated by the contour-plots of B^ e, ° shown on both panels. 
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Fig. 12. — Schematic explanation of the situation observed in the pixels where the extrapo- 
lated magnetic field increases upwards. First it increases reaching a maximum at the height 
"HO" plotted in the lower panel of Fig. 11, then the magnetic field decreases again reaching 
the value of the observed chromospheric magnetic field at the height "HI" 



Fig. 13. — Top: 3D view of the magnetic field of the filament on July the 3 . The extrapo- 
lation results use the photospheric magnetic field as lower boundary condition. The sample 
of field lines is plotted in the filament region to outline its structure. Top background image: 
vertical component of the magnetic field at the bottom of the extrapolation cube. Purple 
and blue: B^ l ' E > 0; green and red: B^ l ' E < 0. Bottom: Extrapolation results for the 3 rd 
of July using the chromospheric magnetic field as lower boundary condition. Bottom back- 
ground image: vertical component of the magnetic field at the bottom of the extrapolation 
cube B^ e,E (same color scale as in the top panel). 



